Decoherence in quantum walks on the hypercube 
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We study a natural notion of decoherence on quantum random walks over the hypercube. We prove that 
this model possesses a decoherence threshold beneath which the essential properties of the hypercubic quantum 
walk, such as linear mixing times, are preserved. Beyond the threshold, we prove that the walks behave like 
their classical counterparts. 
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I. INTRODUCTION 

The notion of a quantum random walk has emerged as an 
important element in the development of efficient quantum al- 
gorithms. In particular, it makes a dramatic appearance in the 
most efficient known algorithm for element distinctness [1]. 
The technique has also provided simple separations between 
quantum and classical query complexity [2], improvements in 
mixing times over classical walks [3, 4], and some interesting 
search algorithms [5, 6]. 

The basic model has two natural variants, the continuous 
model of Childs, et al. [7], on which we will focus, and the 
discrete model introduced by Aharonov, et al. [8]. We refer 
the reader to Szegedy's [9] article for a more detailed discus- 
sion. In the continuous model, a quantum walk on a graph G 
is determined by the time-evolution of the Schrodinger equa- 
tion using kL as the Hamiltonian, where L is the Laplacian of 
the graph and k is a positive scalar to which we refer as the 
"jumping rate" or "energy". In addition to being a physically 
attractive model, it has been successfully applied to some al- 
gorithmic problems as indicated above. 

Such walks have been studied over a variety of graphs 
with special attention given to Cayley graphs, whose algebraic 
structure has provided immediate methods for determining the 
spectral resolution of the linear operators that determine the 
system's dynamics. Once it had been discovered that quan- 
tum random walks can offer improvement over their classical 
counterparts with respect to such basic phenomena as mixing 
and hitting times, it was natural to ask how robust these walks 
are in the face of decoherence, as this would presumably be 
an issue of primary importance for any attempt at implemen- 
tation [10, 11, 12]. 
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In this article, we study the effects of a natural notion of de- 
coherence on the hypercubic quantum walk. Our notion of de- 
coherence corresponds, roughly, to independent measurement 
"accidentally" taking place in each coordinate of the walk at 
a certain rate p. We discover that for values of p beneath a 
threshold depending on the energy of the system, the walk 
retains the basic features of the non-decohering walk; these 
features disappear beyond this threshold, where the behavior 
of the classical walk is recovered. 

Moore and Russell [4] analyzed both the discrete and the 
continuous quantum walk on a hypercube. Kendon and Tre- 
genna [13] performed a numerical analysis of the effect of 
decoherence in the discrete case. In this article, we extend 
the continuous case with the model of decoherence described 
above. In particular, we show that up to a certain rate of de- 
coherence, both linear instantaneous mixing times and linear 
instantaneous hitting times still occur. Beyond the threshold, 
however, the walk behaves like the classical walk on the hy- 
percube, exhibiting ©(nlogn) mixing times. As the rate of 
decoherence grows, mixing is retarded by the quantum Zeno 
effect. 



A. Results 

Consider the continuous quantum walk on the n- 
dimensional hypercube with energy k and decoherence rate 
p, starting from the initial wave function = |0)® n , corre- 
sponding to the corner with Hamming weight zero. We prove 
the following theorems about this walk. 

Theorem 1. When p < 4k, the walk has instantaneous mixing 
times at 

n(2nc — arccos(p 2 /8fc 2 — 1)) 

tmix = , ,, == 

s/l6k 2 -p 2 

for all c G Z, c > 0. At these times, the total variation distance 



2 



between the walk distribution and the uniform distribution is 
zero. 

This result is an extension of the results in [4], and an 
improvement over the classical random walk mixing time of 
0(nlog«). Note that the mixing times decay with p and dis- 
appear altogether when p > 4k. Further, for large p, we will 
see that the walk is retarded by the quantum Zeno effect. 

Theorem 2. When p < 4k, the walk has approximate instan- 
taneous hitting times to the opposite corner (1, ... ,1) at times 

_ 2nn{2c+l) 

for all c G Z, c > 0. However, the probability of measuring an 
exact hit decays exponentially in c; the probability is 



Pint = 



;«r(2r+l) 



2 + 2 



In particular, when no decoherence is present, the walk hits at 
thit = " a ^2^ +1 ^ ' an d it does so exactly, i.e. P/„> = 1. For p > 4k, 
no such hitting occurs. 



B. The continuous quantum walk on the hypercube 

A continuous quantum walk on a graph G begins at a dis- 
tinguished vertex vo of G, the initial wave function of the walk 
being ^q, where (^o| v ) = 1 if v = v and otherwise. The 
walk then evolves according to the Schrodinger equation. In 
our case, the graph is the n-dimensional hypercube. Con- 
cretely, we identify the vertices with n-bit strings, with edges 
connecting those pairs of vertices that differ in exactly one bit. 
Since the hypercube is a regular graph, we can let the Hamilto- 
nian H be the adjacency matrix instead of the Laplacian [14]; 
the dynamics are then given by the unitary operator U, = e' Ht 
and the state of the walk at time t is *P ( = U^q. 

The following analysis makes use of the hypercube's prod- 
uct graph structure; this structure will be useful again later 
when we consider the effects of decoherence. The analysis 
below diverges from that of Moore and Russell [4] only in 
that we allow each qubit to have energy k/n instead of \jn. 
The energy of the entire system is then k. Let 

_ f k/n\ 
aX ~\k/n )' 

and let 



This result is a significant improvement over the exponen- 
tial hitting times of the classical random walk, with the caveat 
that decoherence has a detrimental effect on the accuracy of 
repeated hitting times. 

Finally, we show that under high levels of decoherence, the 
measurement distribution of the walk actually converges to the 
uniform distribution in time ©(« log«), just as in the classical 
case. 

Theorem 3. For a fixed p > 4k, the walk mixes in time 
©(nlogn). 

In the remainder of the introduction, we describe the contin- 
uous quantum walk model, and recall the graph product anal- 
ysis of Moore and Russell [4]. In the second section, we de- 
scribe our model of decoherence, derive a superoperator that 
governs the behavior of the decohering walk, and prove that 
it is decomposable into an n-fold tensor product of a small 
system. We then fully analyze the small system in the third 
section, and use those results to draw conclusions about the 
general walk in 3 distinct regimes: p < 4k, p = 4k, and p > 4k. 
These regimes are roughly analogous to underdamping, criti- 
cal damping, and overdamping (respectively) of a simple har- 
monic oscillator with damping rate p and angular frequency 
2k. 



ff=J> 



where the j'th term in the sum has <J X as the y'th factor in the 
tensor product. Then we have 



U t = e m = J] 1 



1 = [e 1 ^] 



0n 



cos(kt/n) i sin(fo/n) 
i sin(kt/n) cos(kt/n) 



Applying U t to the initial state 'Pq = |0)®", we have 



U,V = 



cos 



kt 



|0)+/ sin( - ) |1 



which corresponds to a uniform state exactly when ^ is an 
odd multiple of j . 

H. A DERIVATION OF THE SUPEROPERATOR 

We begin by recalling a model of decoherence commonly 
used in the discrete model, with the intention of deriving a 
superoperator U t , acting on density matrices, which mimics 
these dynamics in our continuous setting. The discrete model, 
described in [13], couples unitary evolution according to the 
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discrete-time quantum random walk model of Aharonov et 
al. [8] with partial measurement at each step occurring with 
some fixed probability p. Specifically, the evolution of the 
density matrix can be written as 

pt+i = (l-ppptU^+p^PiUptU^Pi 

i 

where U is the unitary operator of the walk, i runs over the 
dimensions where the decoherence occurs, and the P, project 
in the usual "computational" basis [13]. 

In the continuous setting, the unitary operator that governs 
the non-decohering walk is Ut = e~' Ht , where H is the normal- 
ized adjacency matrix of the hypercube times an energy con- 
stant. To extend the above decoherence model to this setting, 
recall that the superoperator U t ® u} associated with these dy- 
namics has the property that 

dUt ® U ' = i (e- im ® e mt ) [1 ® H-H® t] ; 
dt y ' 

wishing to augment these dynamics with measurement occur- 
ring at some prescribed rate p, we desire a superoperator S t 
that satisfies 

St+dt = S, \e- mdt ® e mdt ] [(1 - pdt)t + pdt(P)] 

where P is the operator associated with the decohering mea- 
surement. Intuitively, the unitary evolution of the system is 
punctuated by measurements taking place with rate p, analo- 
gous to the discrete case. 

Letting e~' Hdt = 1 — iHdt, we can expand and simplify: 

S t+dt = S, [e- iHdt ® e mdt ] [(1 - pdt)t + pdt(P)} 

= S,[(1- iHdt) ® (1 + iHdt)] [(l-pdt)t + pdt(P)] 

= S,[l®l+idt(t®H-H®t)-pdtl®t+ pdt(P)] . 

In terms of a differential equation, 

dSt _ St+dt ~ St 
dt dt 

_ S t [i ® 1 + idt{l®H-H® 1) -pdt {I® I +P)] -S t 

dt 

= St[i{t®H-H®i)-pt®l+p(P)) . 

The solution is 

S t =exp([i(l®H-H®l)-pl® i+p(P)]t) . (H.l) 

We now define the decoherence operator P. This opera- 
tor will correspond to choosing a coordinate uniformly at ran- 
dom and measuring it by projecting to the computational basis 
{|0), |1)}. Let IIo and IT be the single qubit projectors onto 
|0) and |1), respectively. We define 

p = - E [n^n^+n^ni] 

n KKn 



where n[,= I® - ■ -®t®Yl ®t®- ■ -®t with the nonidentity 
projector appearing in the zth place. We define Yl\ similarly, so 
that TVj ignores all the qubits except the ith one, and projects 
it onto |;') where 6 {0, 1}. Note that 

TVj ® n^. = [i ® i] ® ■ ■ ■ ® [n^ ®Yij}®---®{n®ii} 

for j G {0,1}. 

A. The superoperator as an n-fold tensor product 

The pure continuous quantum walk on the n-dimensional 
hypercube is easy to analyze, in part, because it is equivalent 
to a system of n non-interacting qubits. We now show that, 
with the model of decoherence described above, each dimen- 
sion still behaves independently. In particular, the superopera- 
tor that dictates the behavior of the walk is decomposable into 
an n-fold tensor product. 

Recall the product formulation of the non-decohering 
Hamiltonian 

u 

H = £ l®---®a x ®---®l 

where 

/ k/n\ 
ax= \k/n ) 

with a x appearing in the y'th place in the tensor product. We 
have given each single qubit energy k/n, resulting in a system 
with energy k. This choice will allow us to precisely describe 
the behavior of the walk in terms of the relationship between 
the energy of the system and the rate of decoherence. 

We can write each of the terms in the exponent of the su- 
peroperator from (II. 1) as follows: 

n 

1® H = £ [1 <g> 1] ® • • • ® [1 ® (7j ® • • • <gl [1 ® 1] , 
7=1 
n 

#®l = £[l®l]®...®[or*<g)l]<g>---<g>[l<gil] . 
Our decoherence operator can also be written in this form: 

p = -£ [nj^+iri^iTi] 

i n 

= - Ed 1 ® i]® •••<£> [n ®no]<g>---<e> 
+[i<g)ii](g>---<g)[ni®ni]®---®[i®i]) . 

The identity operator has a consistent decomposition: 1 ® 1 = 
\ L"=i [1 ® 1] ® • • • ® [1 ® 1]. We can now put these pieces 
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together to form the superoperator: 

S, = exp(i'f (1 <g>#) - U{H® 1) -ptt <g> 1 + p*P) 
= exp f £[1 ® 1] ® • • • ® A® • • • [1 <g> 1] \ 



and thus at time t we have 



=n[i®i] 

7=1 



where 

A = -[(l®i'n<7 r )-(mcr c ® l)-p(t® 1) 
« 

+/j(nig)ni)+p(n ®n )] 

/ ik -ik \ 

t ik -p -ifc 

n —ik — p ik 

\ -iife ifc / 

Notice that for p = 0, [e A ] " = [e- Ua * ® e Ua *]® n , which is 
exactly the superoperator formulation of the dynamics of the 
non-decohering walk. 



III. SMALL-SYSTEM BEHAVIOR AND ANALYSIS OF 
THE WALK 

So far we have shown that the walk with decoherence is 
still equivalent to n non-interacting single-qubit systems. We 
now analyze the behavior of a single-qubit system under the 
superoperator e A . The structure of this single particle walk 
will allow us to then immediately draw conclusions about the 
entire system. 

The eigenvalues of A are 0, — - 



^ and =^L. Here 



a = \J p 2 — 16k 2 is a complex constant that will later turn out 
to be important in determining the behavior of the system as 
a function of the rate of decoherence p and the energy k. The 
matrix exponential of A in this spectral basis can then be com- 
puted by inspection. To see how our superoperator acts on a 
density matrix po, we may change po to the spectral basis, ap- 
ply the diagonal superoperator to yield p,, and finally change 
p t back to the computational basis. At that point we can apply 
the usual projectors Ilo and IT i to determine the probabilities 
of measuring or 1 in terms of time. 

Let *P = 1 0) and p = |¥o) (W 1 . In the diagonal basis, 



Po 



1/2 




a' 
a >. 



e po 



1/2 


, -pt-ai 

i e ^r-(-l + £) 

i -pt+at , „ , 



If we then change back to the computational basis and project 
by Ilo and IIj , we may compute the probabilities of measuring 
and 1 at a particular time t : 



—pt-at -pt+at 

2 + e^ - (1 -p/a)+e~*(I +p/cc] 



-pt-at -pt+at 

2-e^s~(l -p/a)-e^s~(l+p/a) 



which can be simplified somewhat to 



p[o] = - + -*-& 



cos 



p . {fit 

— sin — 

P V2« 



Here we have let j3 = — zee = y 16k 2 — p 2 for simplicity. A 
quick check shows that when p = 0, P[0] = cos 



2(kf 



and 



P[l] = sin 2 (jf), which are exactly the dynamics of the non- 
decohering walk. The probabilities for this non-decohering 
case are shown in Figure 1. 




20 30 40 50 



FIG. 1: The p = case - no decoherence: a plot of P[0] and P[l] 
versus time, for k = 1, n = 5, p = 

The three regimes mentioned before are immediately ap- 
parent. For p < 4-k, j3 is real. When p = 4k, we have j5 = 0, 
which appears to be a serious problem at first glance. Finally, 
for p > 4k, /3 is imaginary. We now address each of these 
three situations in detail. 
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A. The case p <4k : linear mixing and hitting times 

When p < 4k, we recover the perhaps most interesting fea- 
ture of the non-decohering walk: the instantaneous mixing 
time is linear in n. To exactly determine the mixing times for 
our decohering walk, we solve P[Q] = P[l] 
to determining when 



I; this amounts 



1 _gt 
2 



Pt\ P ■ fPt 



equals zero. Clearly the exponential decay term results in mix- 
ing as t — > °°; our principle concern, however, is with the pe- 
riodic mixing times analogous to those of the original walk. 
We thus ignore the exponential term when solving the equal- 
ity y=0, which yields 



El 
P 2 



1 +cos(j3f/«) 
1 — cos()3f/n) ' 



This equation actually has more solutions than the one we 
started with, because of the use of half-angle formulas for sim- 
plification. The solutions that we want are 



2%c — arccos 



P 

8k 2 



where c ranges over the positive integers. Evidently, the mix- 
ing times still occur in linear time; an example is shown in 
Figure 2. Note also that if we let p = 0, we have t m \ x = 
nn(2c — 1)/ (4k), which are exactly the nice periodic mixing 
times of the non-decohering walk. In the decohering case, 
however, these mixing times drift towards infinity, and cease 
to exist altogether beyond the threshold of p = 4k. This proves 
Theorem 1 . 

We now wish to determine when our small system is as 
close as possible to Since our large-system walk begins 
at jO}®", this will correspond to approximate hitting times to 




the opposite corner These times correspond to local 

maxima of Pfll; the solutions are 



thit = 2nn 



2c+l 



where c ranges over the non-negative integers. At these points 
in time, the value of P[l] is 

1 1 -(2c+l)f 

2 2 

which immediately yields Theorem 2. 

B. The breakpoint case p = 4k 

We first observe that t m jx —> °° as p 4k. Hence, we do not 
expect to see any mixing in this case. To analyze the proba- 
bilities exactly, we take the limit of y as p — > 4k. The solution 



lim y 

p^4k 



1 



2kt 
n 



(III.l) 



Indeed, since k, t and n are all positive, y is zero only in the 
limit as t — ► °o. The linear mixing and hitting behavior from 
the previous section has entirely disappeared. As in the critical 
damping of simple harmonic motion, a small decrease in the 
rate p can result in drastically different behavior, in this case 
a return to linear mixing and hitting. We leave the limiting 
mixing analysis of this case for the next section, where we 
develop some relevant tools. 



C. The case p > 4k and the limit to the classical walk 




FTG. 3: The p > 4k case: a plot of P[0] and P[\\ versus time, for 
k= \,n = 5,p = 9 



FIG. 2: The p < 4k case: a plot of P[0] and P[l] versus time, for 
k=\, n = 5,p = 0.5 



The goal of this section is to show two interesting con- 
sequences of the presence of substantial decoherence in the 
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quantum walk on the hypercube. First, we will show that for 
a fixed p > 4k, the walk behaves much like the classical walk 
on the hypercube, mixing in time @(«log«). Second, we show 
that as p — > oo, the walk suffers from the quantum Zeno effect. 
Informally stated, the rate of decoherence is so large that the 
walk is continuously being reset to the initial wave function 
|0)® n by measurement. 

1. Recovering classical behavior 

Consider a single qubit. Let P be the distribution obtained 
by full measurement at time t, and U the uniform distribution: 



And hence, by (III.2), 



P(0) 



1 



where 



7= 4 



Y, PW = \-Y, 17(0) = 17(1) = i 



(-p-a)t (a-/j)r 

! 2n (i - pj a) + e^~ (1 + p/a) 



For x = [x\ , . . . ,x„) £ Zj, 

P n (x) = f[P(xi) and U" (x) = 2~" 

i=l 

are the analogous product distributions in the n-dimensional 
case. To analyze the limiting mixing behavior of the walk, 
we will consider the total variation distance \\P n — U n \\ = 
Y.x\P"(. x ) ~ U n (x)\ between these distributions. In order to 
give bounds for total variation, we will use Hellinger dis- 
tance [15], defined as follows: 



We will make use of the following two properties of Hellinger 
distance: 



1 -H(A'\B") 2 = (1 -H(A,Bf) n 



and 



||A-B|| <2H{A,B) <2\\A-B\\ l l 2 . 



(III.2) 



The first property makes it easy to work with product distribu- 
tions. The second gives a nice relationship between Hellinger 
distance and total variation distance. In our case, 



H(P",U") 2 = 1 - (1 -H(P,U) z ) n 



1- ( -^l+2y+-^Ty 



= i-(i-y+o(r 3 ) 



| 2 <4-4( l-L.+Oif) 



Consider the walk with decoherence rate p > 4k. We have 
a = sj p 2 — 16k 2 < p, where a and p are positive and real. 
It follows that for a fixed p > 4k, y-> and \\P n - U n \\ -»■ 
as t — > oo. Hence the walk does indeed mix eventually, and 
the measurement distribution in fact converges to the uniform 
distribution. Let t = d -«log« where d > is a constant, and 
rewrite / as follows: 



1 



-(p-a) 



— adlogn 

(1 - p/a) + e—r^- (1 + p/a) 



Suppose we choose d such that d > (p — a) . Then y — 
o(n~ l l 2 ), which implies that \\P n - U"\\ = o(l). On the other 
hand, if d < (p — Of) -1 , then y = 0)(n~ 1//2 ) and there exists a 
constant e such that \\P" -U n \\ > e > 0. This shows that the 
walk mixes in time ©(« log«) when p > 4k. Notice that when 
p = 4k, (p — a)^ 1 — (4k) ~ 1 , so that the same technique easily 
extends to that case via equation (III. 1). This completes the 
proof of Theorem 3. 



2. Quantum Zeno effect for large p 

Recall from the previous section that the time required to 
mix when p > 4k is 



t > 



n log 7i 



a 



which clearly increases with p. Further, for large p, p/a tends 
to 1, and hence y tends to 1/2. Notice that y = 1/2 corre- 
sponds to remaining at the initial state forever. We conclude 
that the mixing of the walk is retarded by the quantum Zeno 
effect, where measurement occurs so often that the system 
tends to remain in the initial state. 
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